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The rich variety of iron-based superconductors and their complex electronic structure lead to a wide range 
of possibilities for gap symmetry and pairing components. Here we solve in the two-Fe Brillouin zone the 
full frequency-dependent linearized Eliashberg equations to investigate spin-fluctuations mediated Cooper pair¬ 
ing for LiFeAs . The magnetic excitations are calculated with the random phase approximation on a correlated 
electronic structure obtained with density functional theory and dynamical mean field theory. The interaction be¬ 
tween electrons through Hund’s coupling promotes both the intra-orbital and the inter-orbital magnetic 

susceptibility. As a consequence, the leading pairing channel, conventional s”' , acquires sizeable inter-orbital 

dxy — da; 2 {yz) slnglet pairing with odd parity under glide-plane symmetry. The combination of intra- and inter¬ 
orbital components makes the results consistent with available experiments on the angular dependence of the 
gaps observed on the different Fermi surfaces. 


PACS numbers: 74.20.Pq, 74.70.Xa, 74.20.Rp 

LiFeAs is a stoichiometric superconductor with supercon¬ 
ducting Tc ~ 18 K and no magnetic ordering. [1] De¬ 
spite rather poor nesting [2-5], recent quasiparticle interfer¬ 
ence experiments identify the antiferromagnetic (AF) spin- 
fluctuation mediated mechanism as the predominant pairing 
interaction. [6] ARPES and quasiparticle-scattering interfer¬ 
ence measurements below T], show that the superconducting 
(SC) gaps of LiFeAs are nodeless, with a Fermi surface (FS) 
dependence and a sizable variation along each FS. [2, 7, 8] 
Polarized neutron diffraction as a function of temperature has 
shown a suppression of the local spin susceptibility in the SC 
phase, suggesting singlet pairing. [9, 10] 

In theoretical studies, the AF spin-fluctuation mediated 
pairing [11-14] and a combination of AF spin-fluctuation and 
orbital fluctuation mediated by phonons have been investi¬ 
gated. [15, 16] However, all studies are performed in the one- 
iron unit cell with various unfolding algorithms used to em¬ 
bed the correct symmetry. [17-21] This procedure is exact 
only for computing in-plane pairing. In addition, the SC gap 
equation is usually projected on the FS, the pairing interac¬ 
tion is symmetrized, [11] and the resulting equation is always 
solved in the BCS approximation. All of the above simplifi¬ 
cations must be questioned before we can be confident of the 
results. Furthermore, for Fe-based superconductors (FeSCs) 
with a non-symmorphic point-group, [22] anti-symmetry of 
fermions does not place a constraint on the parity of the SC 
pairing channel. [23, 24] This allows for even-parity — dyz 
inter-orbital pairing [25], or for dxy — dxz{yz) odd parity spin 
singlet pairing when there is orbital weight at the Fermi level 
from orbitals with different in-plane mirror reflection symme¬ 
try [26]. 

Hence, here we revisit spin-fluctuation mediated pairing 
by considering both Fe-Sd and As-4p orbitals in the two-Fe 
unit cell. We solve the linearized Eliashberg equations [27] 
in the two-Ee Brillouin Zone (BZ) to investigate SC pairing 
and gap symmetry. Since there is increasing evidence that 


superconductivity does not emerge as a ES instability [40], 
we work in the orbital representation instead of projecting the 
gap equation on the ESs. Our results show that in the lead¬ 
ing channel, with the conventional symmetry, odd parity 
inter-orbital pairing accompanies the usual intra-orbital pair¬ 
ing and increases with interactions, in particular with Hund’s 
coupling. In contrast to previous studies [8, 11-13] we And 
that this state can reproduce the angular dependence of the 
gap on the electron pockets. 

Electronic structure In LiFeAs, the bandwidth observed 
in ARPES is narrower than in LDA calculations and there are 
experimental evidences of long-lived magnetic moments. [9] 
This indicates the importance of correlations, so we em¬ 
ploy the LDAh-DMET method to obtain the electronic struc¬ 
ture. [41-43] Eig. 1 illustrates the LDAh-DMET partial spec¬ 
tral weight, A;;(k,0), of Fe t 2 g- orbitals d^y and dxz,yz on 
the FSs of LiFeAs. [44] The Fe Cg orbitals dz :2 and dx 2 -y 2 hy¬ 
bridize with As-p orbitals and contribute to the spectral weight 
lying above and below the Fermi level. The FS consists of 
three hole-like and two electron-like sheets around the center 
and corners of the BZ respectively. The two inner hole pock¬ 
ets are predominantly composed from dxz and dyz orbitals. 
The smallest hole pocket crosses the Fermi level only in close 
vicinity to the F point. It hybridizes with the dz 2 orbital near 
Z point and is closed there, while remaining 2D away from 
this point. The middle pocket has moderate kz dispersion. 
The large hole-like Fermi surface originates purely from in¬ 
plane dxy orbitals and therefore is 2D without noticeable kz 
dispersion. The electron pockets are made from an admixture 
of dxy, dxz and dyz orbitals. The electron pockets intersect at 
small kz and their order flips, i.e., the inner pocket at = 0 
is the outer pocket at kz = tt jc. 

Comparison to LDA, [27] shows that in LDAh-DMFT: (a) 
The two inner hole pockets shrink while the outer one ex¬ 
pands. (b) The middle hole pocket also deforms and takes 
on a butterfly shape at small kz- [45] (c) At finite kz, the outer 
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FIG. 1. (Color online) Partial spectral weight, Aii(k, 0), of Fe 
t 2 g- orbitals on the FS in the k^-ky plane with = 0 (left), and 
kz = njc (right) obtained from the LDA+DMFT calculation. Here 
the dxy, dxz, and dyz orbitals are illustrated by green, blue and red 
colors, respectively. The ai pocket crosses the Fermi level only in 
close vicinity to the F point (not visible on this scale). 

hole-pocket acquires some dxz and dyz orbital weight in the 
direction of the A point, (d) The shrinkage of the two inner 
hole pockets leads to larger patches where dxz and dyz orbitals 
mix on these pockets, (e) The electron pockets are moderately 
expanded and they become closer to each other. [27] 

The t 2 g orbitals are the most strongly correlated [43, 45] as 
is apparent from the mass enhancements m */ttilda = 2.0, 
1.85, 3.13 and 2.7 for dz^, dx 2 -y 2 , dxy, and dxz,yz orbitals, 
respectively. The dxy orbital has the strongest mass enhance¬ 
ment and shortest quasi-particle lifetime. 

Effective pairing interaction A SC instability in the sin¬ 
glet channel occurs when the corresponding pairing suscep¬ 
tibility diverges as one lowers temperature. A divergent sus¬ 
ceptibility signals the appearance of a pole in the correspond¬ 
ing reducible complex vertex function, which describes all 
scattering processes of two propagating particles. Using the 
Bethe-Salpeter equation, the condition for an instability is that 
an eignvalue of the matrix —becomes unity. Here 
r*rr,s jg j-jjg irreducible vertex function (effective pairing in¬ 
teraction) in the singlet channel, and is the bare suscepti¬ 
bility in the particle-particle (p-p) channel. [27, 46, 47] 

The density/magnetic fluctuations contribute to the pair¬ 
ing interaction by entering the ladder vertex defined by 

ph = Xph^ + Xph^ 

where pirr.m(d) (jgnote respectively the dressed 

susceptibility and the irreducible vertex function in the mag¬ 
netic (density) channel. [27] These vertices can be calculated 
in the DMFT approximation. [48] However, such a calculation 
is prohibitively difficult for multiorbital systems at the low 
temperatures necessary to study superconductivity, [27] hence 
here we employ the random phase approximation (RPA). [49] 
In RPA, the irreducible vertex function is replaced by a static 
effective vertex which is parametrized by the screened intra¬ 
orbital Hubbard interaction, Ug, and the Hund’s coupling 
Js- [16, 27, 50, 51] The inter-orbital interaction and pair hop¬ 
ping are determined assuming spin-rotational symmetry. Note 
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FIG. 2. (Color online) Several components of the pairing interaction 
of LiFeAs at ksT — 0.01 eV in the particle-hole channel. There 
are two sets of screened interaction parameters yielding the same 
magnetic Stoner factor, namely Js = O.lUs, Us = 2.4 eV on the top 
and Js = O.SUs, Us = 1.68 eV on the bottom. The legend for the 
color coding is spread over both figures. 


that even though the static effective vertices Ug and Jg capture 
Kanamori-Briickner screening effects, they do not fully cap¬ 
ture the dynamics of screening. In particular, the RPA treat¬ 
ment misses the fact that at high fermionic frequencies one 
should recover the bare interactions. 

Fig. 2 shows the pairing interaction, Ilph, at = 

0.01 eV for two sets of screened interaction parameters that 
yield the same magnetic Stoner factor. [52] Here we only 
present the intra-sublattice components because the inter¬ 
sublattice components are relatively small. In what follows, 
we focus on the Fe-1 and Fe-2 (on A and B sublattices re¬ 
spectively) t 2 g orbitals; dxy will be referred as 2 (7) and dxz 
and dyz orbitals as 4 (9) and 5 (10). The dominant effec¬ 
tive pairing interaction components are repulsive. As can be 
seen in Fig. 2(a), due to better nesting, the dxy intra-orbital 
(22; 22) pairing vertex is dominant and the dxz(yz) intra¬ 
orbital (44; 44) is sub-dominant, yet on average it is larger 
than inter-orbital vertices (22; 44) and (44; 55). 

However, at larger Jg/Ug the situation changes. For a 
fixed Stoner factor (proximity to magnetic transition) upon in¬ 
creasing the Jg/Ug ratio from Fig. 2(a) to Fig. 2(b), the dxy 
intra-orbital pairing component decreases while the dxz(yz) 
intra-orbital components and the inter-orbital components in¬ 
crease slightly. This shows that a higher Jg, through cou¬ 
pling to the more correlated dxy orbital, compensates the 
decrease of spin susceptibility expected from the lower Ug 
(Fig. 2(b)). [27] Furthermore, since Hund’s coupling corre¬ 
lates different orbitals, the inter-orbital components increase. 
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FIG. 3. (Color online) Real part of the several intra-sublattice com¬ 
ponents of the generalized particle-particle bare susceptibility at the 
lowest fermionic/bosonic Matsubara frequency. 


becoming comparable with the dxz{yz) intra-orbital compo¬ 
nents. The dxy intra-orbital vertex becomes less dominant at 
larger Jg/Ug- [53] This behavior of the magnetic susceptibil¬ 
ity reflects itself directly in the pairing interaction (see sup¬ 
plemental material for the dressed susceptibilities in magnetic 
and charge channels). 

Bare particle-particle susceptibility The generalized bare 
susceptibility in the p-p channel also enters the gap equa¬ 
tion. [27] Fig. 3 shows the real part of several compo¬ 
nents of the generalized p-p bare susceptibility at the lowest 
fermionic/bosonic frequencies. The intra-orbital components 
are purely real. Both real and imaginary parts (see SM) show 
peaks at the position of FSs. For example, going from the F to 
the X point in the top panel, the three peaks are respectively 
related to the inner hole pocket with dxz weight in close prox¬ 
imity to F, the middle pocket with dyz weight and the outer 
pocket with dxy weight. The peak heights are directly pro¬ 
portional to the corresponding orbital weight on the FSs and 
inversely proportional to the Fermi velocity. The peak widths 
are induced by correlation effects, implying that electrons near 
FSs may contribute to the Cooper pairing. In a non-interacting 
system the peak widths go to zero at zero temperature. [54] 
The larger 22; 22 peak component in the M — T direction, 
compared with the M — X{Y) direction, indicates that the 
SC gap on the outer electron pocket is larger in the M — T 
direction. 

In the BCS approximation, only real parts survive for the 
components considered here, due to a summation over Mat¬ 
subara frequencies. In this case, the inter-orbital pairing is 
suppressed. Including the imaginary part in the full gap equa¬ 
tion changes this trend. The imaginary parts of the inter¬ 


orbital components change sign between corner and center of 
the BZ. They have some symmetries that transfer to the gap 
function: (i) They are odd under exchange of orbital indices, 
(ii) There is also a tt phase difference between the two Fe ions 
(see SM). 

SC pairing symmetry in LDA-\-DMFT+RPA The leading 
pairing channel is a channel with dominant dxy, dxz and 
dyz intra-orbital pairing. In our gauge, the gap function 
components have both real and imaginary part which sat¬ 


isfy 


= -ImA 


AA(BB) 
U 


All intra-orbital com¬ 


ponents change sign between center and comer of the BZ 
(see Fig. 4), as expected in conventional pairing. The 
dxy intra-orbital component dominates, but has a small value 
on the 7 pocket. The dxz and dyz intra-orbital components 
are out of phase, i.e. ~ (not shown). 

They take large values on the a 12 hole pockets. The inter¬ 
sublattice components are much smaller than intra-sublattice 
ones, The largest inter-sublattice 

component is X^^f. In orbital basis, the gap functions do 
not change much between kz — 0 and kz = tt / c , hence we 
present only kz = 0 results. 

In agreement with the above pairing-interaction analy¬ 
sis, upon increasing Jg/Ug the dxz/yz intra-orbital pairing 
strengthen. Furthermore, the dxy-dxz and dxy-dyz inter¬ 
orbital pairings increase. Although they vary on a smaller in¬ 
terval, they are comparable with the dxz/yz intra-orbital com¬ 
ponents on the electron FSs (compare Fig. 4 top and bottom 
panels). 

We verify that the gap function components of the lead¬ 
ing channel satisfy the relations XfJ‘^^^^\'k,iu!m) = 




and 


j,AA(BB),, . 


^AA{BB) 


(—k, —iWm). [55] The first relation says that 
the superconducting state does not break parity: In 
LiFeAs the inversion center is located in the middle 
of Fe-Fe link. Under parity operation the sublattice A 
maps to sublattice B and vice versa and k ——k. The 
components of the gap function also satisfy the rela- 


tio” ^tp^^’^\kx,ky,iu}m) = pi^pi^Xf^f}'^^\k 


X, ky ,: 




where pi denotes the parity of orbital I with respect to in-plane 
mirror reflection symmetry. [56] This symmetry is defined 
by in-plane mirror reflection followed by a half-translation, 
expressed in units of the two-Fe unit cell, {cr^|H0}. Thus, 
the intra-orbital components on the two Fe are equal, while 
the inter-orbital components between one even-parity (dxy) 
and one odd-parity (dxz, dyz) orbital, change sign between 
two Fe-ions. These components are the parity-odd under 
{(T^|^|0} spin singlet pairings. [26] Furthermore, as can 
be seen from Fig. 4, the in-plane intra-orbital components 
satisfy Xfi^^^^\kx,ky) = Xfi^^^^\-kx,-ky), while 

the inter-orbital components between dxy and dxz{yz) 

satisfy Xf^f^^^\kx,ky) = -Xf^f^^^\-kx,ky) or 

Xt^f^^^\kx,ky) = -Xt^f^^^\kx,-ky) . 

Our calculations show that the gap symmetry of the lead¬ 
ing channel is conventional . Indeed, although there is a 
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FIG. 4. (Color online) Top panels: The real part of the d^y (left) and 
dxz (right) in-plane intra-orbital components of the SC gap function 
at the lowest Matsubara frequency with largest eigenvalue in the or¬ 
bital representation for Js/(7s =0.3 and ksT = 0.01 eV. The imag¬ 
inary part can be obtained from ImAj; = —ReA;;. Bottom panels: 
The real/imaginary part of the inter-orbital components of the SC 
gap function on sublattice A in the orbital representation, . The 
corresponding components on sublattice B are out of phase with the 
displayed components, i.e. A®® = — A(J^. The lines show one 
quarter of the Fermi surfaces. 


phase difference between the dxz and dyz components of the 
gap function in the orbital basis, this phase difference is re¬ 
moved by another phase difference that arises when going to 
the Bloch basis corresponding to the ai 2 pockets. [27] In the 
subleading pairing channel, the dxy intra-orbital component 
is in phase with dyz and out of phase with dxz intra-orbital 
components, which in the band representation gives gap 
symmetry with a sign change between 0 : 1.2 and 7 pockets 
and between electron pockets and accidental nodes on the (32 
pocket. [14] 

Finally, we comment on the SC gap magnitude on different 
FSs. [57] Diagonalizing the Bogoliubov quasi-particle Hamil¬ 
tonian leads to a gap magnitude which has predominant cos 40 
angular dependence on all pockets, as can be seen from Fig. 5. 
The angular dependence of the gap on the 7 and of the aver¬ 
age gap on the /3i 2 pockets are consistent with ARPES data: 
The gap is maximum at 0 = 0,7r/2 and decreases when ap¬ 
proaching 0 = 7r/4 (the direction toward M-point) on the 7 
pocket, while the average gap is maximum at 0 = 7r/4 (di¬ 
rection toward F-point) on the (3 pockets and decreases when 
approaching 0 = 0,7r/2 where the two pockets cross. The 
gap on the (32 electron pocket is increased in the direction of 
F-point due to a larger dxy orbital content with a large pair¬ 
ing amplitude (see Fig. 4, upper panels). The gap on the /3i 
electron pocket also shows a local enhancement at 0 = 7r/4. 
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FIG. 5. (Color online) For Jb/Us = 0.3, the SC gap magnitude (in 
units of the average gap magnitude on the ai pocket) as a function 
of the angle 6 measured at the F and M points with respect to the x 
axis for kz = 0 FSs. 


Due to interchange of electron pockets as a function of kz, the 
gap on the inner pocket becomes larger than that on the outer 
pocket at a finite kz- Hence, for these pockets, a direct com¬ 
parison with ARPES data has to take averaging over a range 
of kz into account. [58] The ratio between the average gap 
magnitude on (3 pockets and 7 pocket is also consistent with 
ARPES results [7, 8]. However, the gap magnitude on the a 
pockets is not the largest. This discrepancy with ARPES re¬ 
sults may come from the fact that ARPES is performed at very 
low temperature while the linearized Eliashberg gap equation 
is valid at temperatures infinitesimally close to the transition 
temperature. The tunneling spectroscopy study of LiEeAs has 
shown a temperature evolution of superconductivity. [59] A 
calculation at a lower temperature shows that the sharp peaks 
in the 44 and 55 bare paring susceptibilities, Eig. 3(a), grow 
faster than the wider peak for 22. This leads to an increase of 
the gap on the a pockets at lower temperatures. 

Conclusion Solving the full linearized Eliashberg gap 
equation with both real and imaginary parts and including cor¬ 
relations in the LDAh-DMET framework leads to a detailed de¬ 
scription of the leading pairing channel in LiEeAs. Account¬ 
ing for correlations in the spin fluctuation approach allows to 
correctly capture not only nesting effects but also Fe-d orbital 
fluctuating moments with orbitally dependent dynamics. Al¬ 
though the intra-orbital dxy spin susceptibility is dominant, 
Hund’s coupling between orbitals on individual Ee atoms pro¬ 
motes both the intra-orbital dxz{yz) component and the inter¬ 
orbital dxy-dxz(yz) components of the magnetic susceptibil¬ 
ity. As a consequence, the leading paring channel, conven¬ 
tional 5“' , acquires inter-orbital singlet pairing component 
with odd parity under glide-plane symmetry. This type of pair¬ 
ing may also be realized in other iron-based superconductors. 
Antiphase s"' pairing [14] is sub-leading. The combination 
of inter-orbital odd-parity and intra-orbital even parity singlet 
pairing leads to a description of the angle-dependence and of 
the relative magnitudes of the gap on the (3 and 7 Eermi sur¬ 
faces that is consistent with state of the art experiments. 
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Supplemental Materials: Correlation-Enhanced 
Odd-Parity Interorbital Singlet Pairing in the 
Iron-Pnictide Superconductor LiFeAs 

Here, we present some details on the electronic structure of 
the LiFeAs in the first section. The Eliashberg equation and 
details of the calculation for the pairing interaction in the ran¬ 
dom phase approximation (RPA) and in the dynamical mean 
field theory (DMFT) are presented in the second section. We 
show how this calculation differs from the BCS approxima¬ 
tion to these equations. Our results on the bare and dressed 
susceptibilities are presented in fourth section. Finally, we 
discuss some symmetries of the Green’s function and their 
consequences on the susceptibilities. 


ELECTRONIC STRUCTURE 

LiFeAs crystallizes in a tetragonal strucmre with a space 
group PA/nmm. Fig. SI shows the crystal structure of 
LiFeAs. In our study the crystal structure is fixed to the ex¬ 
perimental structure. 

We preformed a fully self-consistent LDAh-DMFT calcula¬ 
tion. [SI] We use a local Coulomb integral U — Fq = 4.0 eV 
and a Hund’s coupling J = {F 2 + F4)/14 = 0.8 eV, where 
Ffc are Slater integrals. [S2] We take the fully localized limit 
double-counting correction. We use 32 x 32 x 16 k-points and 
a rotationally invariant interaction matrix. 

The momentum-resolved spectral function, TrA(k, w), 
(trace over all orbitals) along high-symmetry directions is dis¬ 
played in Fig. S2, and compared with the LDA band stmcture. 
Compared to other pnictide families, LiFeAs has much shal¬ 
lower hole pockets around F. The flat top of these pockets im¬ 
plies a large density of states, [S3] which has been proposed 
to promote ferromagnetic fluctuations and triplet p-wave pair¬ 
ing. [S4] However, experiment supports singlet-pairing. 

As can be seen from TrA(k, to), correlation effects lead to 
bandwidth reduction for the Fe-3d states near the Fermi level: 
both the hole-pocket near the F point and the electron-pocket 
near the M-point are pushed towards the Fermi energy. The 
outer hole pocket becomes more incoherent. 

The two low-energy hole-like bands at the F point form a 
degenerate doublet with the two-dimensional irreducible rep¬ 
resentation Eg. [S5] The subscript g and u denote the spatial 
even and odd parity respectively. The two degenerate eigen¬ 
states are predominantly given by (ici — C 2 {dy^ + 

dyj,ic2{d^z + d,^z) Fci{dy^ Fdy^)) where A, B refer toFe 
ions and ci, C 2 are real coefficients. Away from the F point, 
the symmetry reduces and the doublet splits. Depending on 
the k-point, the eigenstate is composed predominantly from 
dxz, or dyz orbital, but a phase factor i between expansion 
coefficients is present everywhere. Along the high-symmetry 
lines T-M and F-A, the Eg symmetry reduces to one A 2 and 
one B 2 state. 

The outer hole-like band at the F-point instead has Big 
symmetry with the d^y + d^y Fe states. Along the high- 
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FIG. SI. (Color online) Crystal structure of LiFeAs. The structure is 
tetragonal with space group PA/nmm. Iron atoms are represented 
by red spheres, arsenic by yellow spheres. Light grey balls are larger 
cations, Li, between the layers. The c-axis points up. Note that a 2 x 
2x2 super-cell is plotted here. The FeAs layers are the electronically 
active part of the compound. 
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FIG. S2. (Color online) LDA electronic band-structure (top) and 
LDA-fDMFT momentum-resolved spectral function (bottom). [SI] 
The colors indicate Fe-d and As-p orbital contents: d^y, are 

denoted by red and green respectively, while all other orbitals are 
displayed in blue. 




FIG. S3. (Color online) Fe t 2 g character of the FSs at k^-ky plan 
with kz = 0 (top panels), and k ^ = tt/c (bottom panels) obtained 
from LDA calculation. Left, middle and right columns show d^y, 
dxz and dyz orbitals respectively. 


symmetry lines T-M and F-X, the symmetry reduces to Bi 
and Ai, respectively. [S5] All the irreducible representations 
at the BZ boundary are two-dimensional [S5] which implies 
that all the Bloch states at any k on the boundary, including 
M-X line, are doubly degenerate. This can be seen from 
Fig. S2. At the M-point, the low-energy doubly degenerate 
electron-like bands are composed of + idy^,dyz — id^^) 
with symmetry A 2 0 B 2 along M-F and E along M-X. The 
doubly degenerate bands with higher binding energy are com¬ 
posed of (d^y + d^y,d^y — d^y) wi^ symmetry Bi 0 B 2 and 
E along M-F and M-X, respectively. [S5] 

Fig. S3 and Fig. S4 illustrate respectively the LDA and 
LDAh-DMFT partial spectral weight, A/i(k, w), of Fe t 2 g or¬ 
bitals on the Fermi surfaces. The two inner hole pockets 
are predominantly composed from dxz and dyz orbitals. The 
smallest hole pocket has strong kz dispersion and is present 
only at small kz, while the middle one has weak kz disper¬ 
sion. The large hole-like Fermi surface originates purely from 
in-plane dxy orbitals and therefore is 2D without noticeable 
kz dispersion. The electron pockets are made from an admix¬ 
ture of dxy, dxz and dyz orbitals with weak kz dispersion. The 
electron pockets intersect at hnite kz and their order flips, i.e., 
the inner pocket at = 0 is the outer pocket at kz = tt/c. 
The orbital character of portions of the (52 pocket oriented to¬ 
wards F-points is dxy, while for those portions oriented to¬ 
ward X-point it is dxz or dyz- The (5i pocket originates made 
predominantly from dxz and dyz orbitals. 

Including interactions with LDAh-DMFT, the two inner 
hole pockets shrink while the outer one expands. The mid¬ 
dle hole pocket also deforms and takes on a butterfly shape at 
small kz- The electron pockets are moderately expanded and 
they become closer to each other. 
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FIG. S4. (Color online) Fe t 2 g character of the FSs in the kj;-ky plane 
with kz = Q (top panels), and k ^ = tt/c (bottom panels) obtained 
from the LDA+DMFT calculation. Left, middle and right columns 
show dxy, dxz and dyz orbitals. The inner hole pocket crosses the 
Fermi level only in close vicinity to the F point and is hard to resolve. 


SUPERCONDUCTING INSTABILITY 

A superconducting instability in the singlet channel occurs 
when the pairing susceptibility in this channel diverges as one 
lowers temperature. A divergent susceptibility signals the ap¬ 
pearance of a pole in the corresponding reducible complex 
vertex function, which describes all scattering processes of 
two propagating particles. Using the Bethe-Salpeter equation, 
the reducible vertex function may be written as a function of 
the irreducible one and of the susceptibility that does not in¬ 
clude vertex corrections. We will call this susceptibility bare, 
even though the Green functions can be dressed. One can then 
check that the condition for an instability is that an eignvalue 
of the matrix —becomes unity. Here is the 

irreducible vertex function in the singlet channel, and is 
the bare susceptibility in the particle-particle channel. [S6, S7] 
This yields a non-Hermitian eigenvalue problem, equivalent 
the linearized Eliashberg equation, defined as 




N 

[xlp{Q = 0)] 1,, 


I3I4 


K'd3U;K"Ml6 




(SI) 


where A/j^ (k, iWm) is the gap function which evolves 
smoothly into the off-diagonal (anomalous) self-energy bel¬ 
low Tc- [S7] Ii,...,Iq are combined ion (or sublattice) 
and orbital indices and we have defined K = {'k,ujrn) as 
momentum-frequency four-vectors. Even though the eigen¬ 


value problem is non-Hermitian, the leading eigenvalues are 
real. [S7] The superconducting transition temperature, can 
be obtained as the temperature when the maximum (dimen¬ 
sionless) eigenvalue of Eq. (SI) becomes unity. The eigen¬ 
value problem in Eq. (SI) involves diagonalization of a matrix 
of size N X -F 2) x where N in the number of k 

points, in the maximum number of positive Matsubara 

frequencies and Norb denotes orbital numbers. 

The generalized bare susceptibility in the particle-particle 
channel is 

(S2) 

where G denotes the fully interacting propagator. A factor 
of I in Eq. (S2) arises due to indistinguishability and should 
be considered to avoid double-counting. [S6] This quantity 
is plotted in Fig. (3) in the main text without the prefactor 
N/ksT. 

One way to see that the eigenvector of Eq. (SI) is re¬ 
lated to the anomalous self-energy is as follows. Con¬ 
sider the one-band case in the BCS approximation. In this 
approximation one can neglect the frequency dependence 
of the gap function and of the pairing interaction, which 
allows to perform the summation on Matsubara frequen¬ 
cies. Using bare propagators in Eq. (S2), one can see that 

l)/4ek = tanh(ei(;/2fcBT)/4ek. Here, np denotes the Fermi 
distribution function and we have used the relation nF(ek) = 
(1 — tanh(ek/2/cBT))/2. Inserting this into Eq. (SI), one 
finds 


tanh(ek'/2A:B7’) . 

-ZAb-' 

2ek' 

= A(r)Ak, (S3) 

which, for A = 1, is the Unearlized BCS gap equation. The 
BCS gap equation has the same form as Eq. (S3) but with Ck 
replaced by the quasi-particle excitation E'k = \/e^ + | Akp. 
Considering the Bogoliubov Hamiltonian, one can see that 
Ak is the anomalous self-energy. 



PAIRING INTERACTION 


Evaluation of the effective pairing interaction is a 

quantum many-body problem. [S6, S7] One way is to employ 
Parquet equations to formally rewrite it in terms of the fully 
irreducible vertex function in the singlet channel. A®®'®’’®, and 
of the vertex ladders in the particle-hole density and magnetic 
channels, as [S8, S46] 
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[^‘^^(K' + K + Q)] 


-K'hU--Khl2 


+ K + Q)].K'wu-,-Ki,i,. (S4) 


where a negative sign difference between terms including 
and the corresponding terms in Ref. [S 6 ] is compen¬ 
sated by a negative sign difference in definition of susceptibil¬ 
ities (see Eq. (S5) below). 

The momentum-frequency appearing as an argument 
shows the transferred momentum-frequency, while the two 
momentum-frequency appearing as indices indicate the 
momentum-frequency of one of the incoming external 
legs and one of the outgoing ones. In other words, 
{Q)]Kiii 2 -,K'i 3 i 4 describes the particle-particle scatter¬ 
ing of electrons in orbitals li, I 2 with momenta/frequencies 
{K + Q,—K) to electrons in orbitals ( 3 , (4 with mo¬ 
menta/frequencies {K' + Q,—K'). Here, we are interested 
in the pairing of electrons at zero center of mass momentum 
and zero Matsubara frequency, i.e., Q = 0 component. 

The vertex {Q)]Khi 2 -,K'i 3 i 4 is irreducible in all two- 

particle channels. It is not singular. Its leading order is given 
by a linear combination of the irreducible density and mag¬ 
netic vertices. The vertex ladder functions that appear in 
Eq. (S4) are important for pairing and are defined as 

K'h^iQ) = (S5) 

where the right-hand side is a matrix multiplication in K 
and orbital indices combined. Each matrix multiplication 
has an implied normalization factor {ksT/N). The func¬ 
tions and Xpf^ respectively, irreducible vertex 

functions and generalized dressed susceptibilities in the den¬ 
sity/magnetic channels. Using the Bethe-Salpeter equation in 
the particle-hole density/magnetic channels, can be de¬ 
composed into bare susceptibility and vertex correction with 
the help of the irreducible vertex functions: 

xfr\Q) = xlh{Q) 

- { + )xfr'^ (Q)xO JQ). (S6) 

Again the second term on the right-hand side is a matrix multi¬ 
plication with the corresponding normalization factor and the 
bare generalized susceptibility is given by 

[Xph{Q)]Khl2-,K'hl4 = -{-j^^)GK+Q,hl3^K,lil2^K,K' ■ 

(S7) 

Eor further progress, one has to do some approximation. 
One can calculate the ladder functions in the particle-hole 
channel, by simply replacing all irreducible vertex 

functions with a static and momentum-independent interac¬ 
tion. This corresponds to the random-phase approximation 


(RPA). In RPA another simplification appears due to time- and 
spatial-locality of the vertex functions. This allows to per¬ 
form the summations over momentum-frequency of the ex¬ 
ternal legs. This greatly simplifies the equations and it also 
reduces the size of matrices, which now are only functions 
of transferred momentum-frequency vector and of orbital in¬ 
dices. Eor example, for the bare susceptibility in the particle- 
hole channel we obtain 






,kj^ 2 
^ N ’ 


[Xph{Q)\Khl2-,K'l3l4 

KK' 


-(^-) y] GK+Q,hl2GKMl2 

K 


(S 8 ) 


Thus, one can drop the external momentum-frequency in the 
Parquet equations and rewrite the pairing interaction at Q = 0 
as 




-- K)]i2W,i3h + 

- + li^phiK' + 

(S9) 


In the DMET framework, the irreducible vertices are still 
local but they become frequency dependent [S9, SIO] and 
in the context of iron pnictide superconductivity have only 
been calculated in the LDAh-DMET study of Refs. [SI 1, S12]. 
However, obtaining the full frequency dependence of the 
DMET vertex is numerically very expensive and challeng¬ 
ing for the multiorbital systems at the low temperatures. A 
calculation of the dressed susceptibility in DMET is possi¬ 
ble by quantum Monte Carlo (QMC) in the Hirsch-Eye algo¬ 
rithm [SI3] or hybridization expansion algorithm [S14]. The 
irreducible vertex can be obtained from the dressed suscep¬ 
tibility using Bethe-Salpeter equation. However, both meth¬ 
ods have serious limitations. The Hirsch-Eye algorithm is 
only applicable at very high temperature, which is not suit¬ 
able for studying superconductivity. In the temperature range 
accessible by this method, the pairing instability is very weak. 
The hybridization expansion algorithm gives good results for 
the dressed susceptibility only at low frequencies. Indeed, 
the noise grows very fast for moderate and high frequen¬ 
cies. [S15] To extract the irreducible vertex from the dressed 
susceptibility, using the Bethe-Salpeter equation, requires all 
frequencies, consequently, even the low frequency data of the 
full vertex are not usable. Eurthermore, even if we assume that 
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we have the frequency dependent irreducible vertex, the size 
of eigenvalue problem in the Eliashberg equation for multi¬ 
band systems is so large that only a limited number of low¬ 
est Matsubara frequencies can be considered. On this lim¬ 
ited number of frequencies, the irreducible vertex does not 
reach its asymptotic behavior. Therefore an RPA analysis 
with static screened irreducible vertices (representing an aver¬ 
age over frequency) seems a better approximation. Hence we 
take local frequency-independent irreducible vertices. [S16] 
Since the pnictides are weakly to moderately correlated sys¬ 
tems, including the DMFT propagators in the resulting RPA- 
like equations gives a better description of Fermi surface nest¬ 
ing, and accounts for the interaction of propagating particles 
with their environment. We think the results should be more 
reliable than standard LDAh-RPA. 

Even with the above approximations for the equa¬ 

tion for r*''’’’® depends on the incoming and outcoming 
momentum-frequency vectors. Thus, full matrices appear in 
the Eliashberg equation. 

To our knowledge, in all the previous studies another ap¬ 
proximation, the so called BCS approximation, has been used 
to solve the Eliashberg equation. In this approximation, since 
the pairing interaction falls off rapidly on a frequency scale 
which is small compared with the bandwidth, one can ne¬ 
glect the frequency dependence of the gap function. Further¬ 
more, in this approximation one takes into account only the 
contribution from the lowest Matsubara frequency, the effec¬ 
tive interaction is symmetrized, and its imaginary part is ne¬ 
glected, [S17] making the Eliashberg equation a real symmet¬ 
ric equation. Then the summation over Matsubara frequencies 
in the bare p-p susceptibility is performed to further simplify 
the Eliashberg equation. Moreover, often the resulting equa¬ 
tion is projected on the Fermi surfaces and is solved assuming 
areal gap function with a specific symmetry. [S8, SI 8] Strong 
coupling treatments such as those of Refs. [S19] and [S20] do 
not restrict the pairing function to the vicinity of the Fermi sur¬ 
face but do not take into account the retardation of the pairing 
interaction. 


Here, we solve the full equation in the physical 2-Fe Bril- 
louin zone with 10-d orbitals at kgT — 0.01 eV. We use a 
16 X 16 X 2 k-mesh and 12 Matsubara frequencies. The con¬ 
vergence of the results have been checked by changing the 
above values to a 8 x 8 x 2 k-mesh and 24 Matsubara fre¬ 
quencies. The three largest eigenvalues of the former set are 
~ 0.125, 0.099, 0.085 and they change by less than 10% for 
the second set. More importantly, the gap symmetries do not 
change between the two sets. The leading channel is conven¬ 
tional s"' , the next-leading channel is anti-phase while 

the third one has dxy symmetry. 


The irreducible vertex function in density/magnetic chan¬ 
nels are defined as: -f 

In RPA, the irreducible vertex function is replaced 
by the antisymmetrized static Coulomb vertex, 
which, in terms of the interacting part of Hamiltonian 

( 1 / 2 ) Sc 


T(T(T' 
t/ J-) 


de¬ 


fined by ~ 

1 ! , where a = —a. For a local interaction the follow- 
ing forms for density/magnetic irreducible vertex functions 
are obtained 


, 0 ,cr (7 


j^irr,d(ni) 
/i ^2 ;^3^4 


Us{Us) 

-C/'+ 2J,([/') 
2 [/' - Js{Js) 

J'siJ's) 

0 


h — h — h — h 

/l = (3 ^ (2 — ^4 
/l = (2 7^ ^3 ~ ^4 

/l = (4 ^ (2 — ^3 
Otherwise 


(SlO) 


where Us and (7' denote the screened static local intra- and 
inter-orbital density-density interactions while Jg and J' are 
Hund’s coupling and pair-hopping interactions. Due to local¬ 
ity of the interaction the four orbital indices belong to same 
ion. We also assume spin rotational invariance so the equal¬ 
ities U'g = Us — \Js [S21] and J' = Jg are satisfied. Fi¬ 
nally, the fully irreducible vertex {Q)\ki 412 -,k' 1314 , is 

replaced by \ transformed to the particle- 

particle channel. 

In our study we used two sets of screened interaction pa¬ 
rameters yielding the same magnetic Stoner factor, namely 
Jg = O.lUs, Us = 2.4 eV and Jg = 0.3C/g, Us = 1.68 eV. 
The values for Us and Jg are relatively standard in the litera¬ 
ture that uses the RPA approach for the pairing vertex. They 
are smaller than cRPA values [S22] in order to account for 
the RPA overestimation of spin fluctuations. We found that 
for a given Jg/t/g ratio changing them within a limited range 
does not change the qualitative aspects of our results for the 
superconducting state. 


RPA SUSCEPTIBILITIES 

The top panels of Fig. S5 compare important components 
of the FDA bare susceptibilities in the particle-hole channel, 
Un = 0) at keT = 0.01 eV, with their FDAh-DMFT 
counterpart along a high-symmetry path. The Fe-1 (Fe-2) <7^2, 
dxy, dx'i-y'i, dxz, dyz otbitals are labeled respectively as 1 (6), 
2 (7), 3 (8), 4 (9), and 5 (10). Due to better nesting, [S23] 
the dominant component of the FDA bare susceptibility is 
the dxy intra-orbital component, [xlhh 2-,22 (= [Xph] 77-,77 = 
[xU 27 ; 27 ), with incommensurate peaks around the M (A) 
point, which come from nesting between the hole and electron 
pockets. The inter-orbital dxy-dxz{yz) component, (24; 24) 
and (25; 25) respectively, is rather large and shows a peak 
centered at the M (A) point. The intra-orbital dxz{yz) com¬ 
ponents (44; 44), and the inter-orbital dxz-dyz components 
(45; 45), have almost the same magnitude in the entire k- 
space. Nesting between the hole (electron) pockets is weak 
and gives small peaks at small momenta. 

Upon introducing correlations within the FDAh-DMFT 
framework, the absolute values of the bare susceptibilities are 
suppressed for all components. However, the AF fluctuations 
remain the leading instability. The small peaks due to intra¬ 
hole (electron) pockets scattering are shifted towards smaller 
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FIG. S5. (Color online) comparison between LDA (panel (a)) 
and DMFT (panel (b)) components of the bare susceptibility, 
[Xph(qi^'" = of LiFeAs at fesT = 0.01 eV in the 

particle-hole channel (top panels). For K summation in Eq. (S8) we 
have used a 32 x 32 x 16 k-mesh and 1024 positive frequencies. 


momenta and smoothed, resulting in a relative increase of zero 
momentum scattering. 

Fig. S6 show several components of the RPA susceptibility 
in the density and magnetic channels for two sets of screened 
interaction parameters that yield the same magnetic Stoner 
factor. The RPA magnetic (charge)-susceptibility is similar 
to the bare one but enhanced (reduced) due to the effect of 
the Coulomb interaction, = Xph/i'^ “ 

The form of the irreducible vertex functions, Eq. (SIO), shows 
that a non-zero Hund’s coupling breaks spin-orbital symme¬ 
try and suppresses orbital fluctuations. This can be seen from 
a suppression of the charge fluctuation upon increasing Jg ( 
compare panels (c) and (d) of the Fig. S6). The magnetic sus¬ 
ceptibility is also clearly larger than the density susceptibility. 
This translates into a repulsive effective pairing interaction. 
This is more pronounced for the intra-orbital components in 
general and in particular for dxy The inter-orbital compo¬ 
nents with for example ( 1,(2 given by (24) 

dxy-dxz or (45) dxz-dyz in Fig. S6, are susceptibilities which 
in RPA are enhanced by inter-orbital Coulomb interaction, as 
can be seen by comparing with Fig. S5(b). 

Interestingly, upon increasing the Js/Ug ratio from 
Fig. S6(a) to Fig. S6(b), the dxy intra-orbital component 
[x^hh 2-22 decreases while the dxy-dxz{yz) and dxz-dyz inter¬ 
orbital components [x)( 5 j] 22 ; 44 ( 55 ) and increase 

slightly. Recall that the calculation is done for a given Stoner 
factor (proximity to magnetic transition). For Jg/Ug = 0.3 
the value of Ug is smaller than for Jg/Ug = 0.1. One would 
have expected that with a smaller Ug intra-orbital magnetic 
correlations would have decreased for all orbitals. Instead, 
only dxy seems strongly affected. This suggests that the larger 
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FIG. S6. (Color online) Several components of the dressed sus¬ 
ceptibility in, respectively, the magnetic channel, density channel of 
LiFeAs at ksT = 0.01 eV in the particle-hole channel. There are 
two sets of screened interaction parameters yielding the same mag¬ 
netic Stoner factor, namely Js = O.lUa, Us ~ 2.4 eV on the left 
and Js = O.SUs, Us = 1.68 eV on the right. Here, Js denotes 
Hund’s coupling and Us the local intra-orbital Hubbard interaction. 
The inter-orbital interaction and pair hopping are determined assum¬ 
ing spin-rotational symmetry. The same line color convention is used 
in all panels. The convention can be read from the insets of parts 
(c),(d) of the figure. 2,4,5 stand respectively for dxy,dxz and dyz. 


value of Jg not only increases inter-orbital spin susceptibility, 
as expected, it also increases intra-orbital spin susceptibility in 
the less strongly correlated orbitals. It is as if, through Jg, spin 
fluctuations in the more correlated dxy orbital increase those 
in the other two t 2 g orbitals. This behavior of the magnetic 
susceptibility reflects itself directly in the pairing interaction. 
Furthermore, due to the negative sign of the corresponding 
off-diagonal dressed susceptibility in the charge channel, the 
magnetic and charge channels for the inter-orbital components 
cooperate to boost the corresponding pairing vertex. 

Bare particle-particle susceptibility The generalized bare 
susceptibility in the p-p channel also enters the gap equa¬ 
tion. Fig. S7 shows the real part in (a, c) and the imaginary 
part in (b, d) for several components of the generalized p-p 
bare susceptibility at the lowest fermionic/bosonic frequen¬ 
cies. The intra-orbital components are purely real. Both real 
and imaginary parts show relatively sharp peaks at the posi¬ 
tion of FSs. Panels (c) and (d) show the inter-sublattice com¬ 
ponents. In the BCS approximation, only real parts survive 
for the components considered here, due to a summation over 
Matsubara frequencies. In this case, the inter-orbital pairing 
is suppressed. Including the imaginary part in the full gap 
equation changes this trend. The imaginary parts of the inter¬ 
orbital components change sign between corner and center of 
the BZ. They have some symmetries that transfer to the gap 
function; (i) They are odd under exchange of orbital indices 
(see Eq. (SI 8)), as can be seen comparing the green and black 
lines or the grey and purple lines in Fig. S7(b); (ii) There is 
also a TT phase difference between the two Fe ions (not shown). 
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FIG. S7. (Color online) Real part (top) and imaginary part (bottom) 
of the several intra-sublattice (left) and inter-sublattice (right) com¬ 
ponents of the generalized particle-particle bare susceptibility at the 
lowest fermionic/bosonic Matsubara frequency. 



FIG. S8. (Color online) For Js/Us = 0.3, the SC gap magnitude (in 
units of the average gap magnitude on the ai pocket) as a function 
of the angle 9 measured at the F and M points with respect to the 
kj; axis for = 0 FSs. Down filled pink triangles are for the P 2 
pocket and up filled blue triangles are for the /3i pocket. The open 
symbols show the gap magnitude at O/tt ~ 0.0,0.25,0.5,, calcu¬ 
lated by setting the inter-orbital components to zero. This calculation 
qualitatively shows that the inter-orbital components reduce the gap 
magnitude where the electron pockets intersect. 


INTER-ORBITAL PAIRING COMPONENTS AND THEIR 
EEFECTS ON THE SC GAP MAGNITUDE 


A quantitative evaluation of the inter-orbital pairing com¬ 
ponent effects on the SC gap magnitude requires a full disen¬ 
tanglement of the changes caused by the inter-orbital compo¬ 
nents. This is not a straightforward task, because inter-orbital 
components are coupled with intra-orbital components in the 
Eliashberg equations. One way to qualitatively understand 
their effects is to remove this component from the gap func¬ 
tion and recalculate the gap magnitude. As can be seen from 
Fig. S8, such an analysis shows that the inter-orbital compo¬ 
nents reduce the gap magnitude, in particular where the elec¬ 
tron pockets intersect, which is where these components are 
maximum. It is worth mentioning that due to anisotropic spin 
fluctuations in our gauge, [S24] the gap function shows a very 
small anisotropy. The presented data here are symmetrized. 


SYMMETRIES OF GREEN’S FUNCTIONS AND 
SUSCEPTIBILITIES 

The normal Green’s function (k, r) describing the 

propagation of electrons from li to (2 with momentum k is 
defined as 


potential. Then complex conjugation yields 
Gr.,,(k,T)- = 

= -(JVcki,.(r)cL,,(0)) = Gi;,,(k,r). (Sll) 
Transforming to Mastubara frequency space we obtain 


(S12) 


This property of the Green’s function implies for the bare sus¬ 
ceptibility in the p-h channel, 

= [x°h(q,-Wn)]cC;ii/2- (S13) 


It also gives a similar relation for the generalized 
bare susceptibility in the p-p channel with zero trans¬ 
ferred momentum/frequency, [App(O)]^ ZaU = 


{k,iuj^)* = [App(O)], , , ik,-iuJm). 


(S14) 

If the time-reversal symmetry is respected then the Green’s 
function transforms as 


r) — (^rCkiiCT('p)c|(.;2^(0)), 

where Ckii<T(T) = exp[(iT-/rA^)T]ckQii<T exp[-(iJ-^Ar)T] 
with r the imaginary time, N = Eka 4zi<TCkZi<T the 
number operator, n the chemical potential, and Tr the time¬ 
ordering operator. Define K = H — fiN, and as the grand 


rGr,;,(k,r) = 

-0(r)e^^Tr[c_kZ3^e-^^cLkz,,,e^-e-^^] 

+ 0(-r)e^^Tr 

= -(T^C_kZ3a(T)cLkZiff(0)) = G1^i^{-k,T) 

(S15) 
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where we have used T| t) = —\.) and T| i) = 
t) with = i. Note that d-orbitals in cubic sym¬ 
metry are time-reversal invariant. It can be seen from their 
definition in terms of Y;™. 

Transforming to Mastubara frequency, the above equality is 
expressed in frequency space as 

~ (S16) 

Combining with the earlier results obtained for complex 
conjugation, we obtain the following identity 

(517) 

which for the generalized bare susceptibility in the p-p chan¬ 
nel, Eq. (S2), gives (assuming SU(2) symmetry) 

= l^^GKMhG-KMh^K,K' 

(518) 
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